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1 Introduction 



There is a general recognition today that basic features of asteroseismology of non-convective final 
stage (FS) stars, such as white dwarfs, neutron stars and and strange quark star, can be properly 
understood working from the model of vibrating solid star, rather than the liquid star model lying 
at the base of asteroseismology of convective main-sequence (MS) stars. In accord with this, most 
of current investigations of an FS-star vibrations rests on principles of solid-mechanical theory of 
continuous media, contrary to the study of the MS-star vibrations which are treated in terms of fluid- 
mechanical theory. This means that super dense matter of FS-stars (whose gravitational pressure 
is counterbalanced by degeneracy pressure of constituting Fermi-matter), possesses elasticity and 
viscosity generic to solid state of condensed matter, whereas a fairly dilute matter of the MS (whose 
internal pressure of self-gravity is opposed by radiative pressure) possesses property of fluidity which 
is generic to the liquid state of a highly conducting condensed matter. This feature of the MS-star 
matter plays crucial role in generation of their magnetic fields in the dynamo processes involving 
macroscopic flows which are supported by energy of nuclear reactions in the central, reactive zone, 
of these stars. In the meantime, in the FS stars, like white dwarfs and neutron stars, there are no 
nuclear energy sources to support convection. The prevailing today view, therefore, is that a highly 
stable to spontaneous decay dipolar magnetic fields of neutron stars are fossil. The extremely large 
intensity of magnetic fields of degenerate solid stars is attributed to amplification of fossil magnetic 
field in the magnetic-flux-conserving process of core-collapse supernova. 

Even still before discovery of neutron stars, it has been realized that, for absence of nuclear 
sources of energy, the radiative activity of these of compact objects should be powered by energy of 
stored in either rotation or vibrations and that the key role in maintaining the neutron star radiation 
should play an ultra strong magnetic field. As is commonly know today, the neutron star capability 
of accommodating such a field is central to understanding pulsating character of magneto-dipole 
radiation of radio pulsars whose radiative power is provided, as is commonly believed, by the energy 
of rigid-body rotation. The discovery of soft gamma-ray repeaters and their identification with 
magnetars [lj - quaking neutron stars endowed with ultra strong magnetic fields experiencing decay 
- has stimulated enhanced interest in the study of models of quake-induced magneto-mechanical 
seismic vibrations of neutron star and resultant electromagnetic radiation. Of particular interest 
in this study are torsional magneto-mechanical vibrations about axis of magnetic dipole moment 
of the star driven by forces of magnetic-field-dependent stresses in a perfectly conducting matter 
and in a permanently magnetized non-conducting matter [2J. Most, if not all, reported up to now 
computations of frequency spectra of poloidal and toroidal Alfven vibration modes in pulsars and 
magnetars rest on tacitly adopted assumption about constant-in-time magnetic field in which a 
perfectly conducting neutron star matter undergoes Lorentz- force-driven oscillations [31 |H [5j El |7J 
US QUI HH H2] (see, also, references therein). A special place in the study of the above Alfven 
modes of pure shear magneto-mechanical seismic vibrations (a-modes) occupies a homogeneous 
model of a solid star with the uniform density p and frozen-in poloidal static magnetic field of both 
homogeneous and inhomogeneous internal and dipolar external configuration, and we start section 
2 with a brief outline of this model. In section 3, a mathematical background of quaking neutron 
star model is briefly outlined with emphasis on the loss of vibration energy caused by depletion of 
internal magnetic field pressure and resulting vibration-energy powered magneto-dipole radiation. 
The decreasing of magnetic field pressure in the star is presumed to be caused by coupling between 
vibrating star and outgoing material which is expelled by quake, but mechanisms of star-envelope 
interaction resulting in the decay of magnetic field, during the time of vibrational relaxation, are 
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Figure V. The lines of magnetic field in fiducial model of neutron star with dipolar external 
and homogeneous internal magnetic field. 

not considered. It is shown that physically meaningful inferences regarding radiative activity of 
quaking neutron star can be made even when detailed mechanisms of depletion of magnetic field 
pressure in the process of vibrations triggered by quakes are not exactly known. This statement is 
demonstrated by a set of representative examples of magnetic field decay. The basic results of the 
model of vibration powered neutron star are summarized in section 4 with emphasis of its relevance 
to astrophysics of magnetars. 



2 Lorentz-force-driven torsion vibrations of a neutron 
star 

To gain better understanding of the basic physics behind the interconnection between seismic and 
radiative activity of quaking neutron star, we start with a brief outline a fiducial model of a solid 
star with frozen-in homogeneous internal magnetic field 

B(r) = Bb(r), B = constant, (1) 
b(r) = [b r = cos 9, be = — sin 9, b^ = 0] 

and dipolar external magnetic field, as is shown in Fig.l. In the last equation, B is the field intensity 
[in Gauss] and b(r) stands for the dimensionless vector-function of spatial distribution of the field. 

The above form of B(r) has been utilized in recent works[ini [HI H2] in which the discrete 
spectra of frequencies of node-free torsional Alfven oscillations has been computed in analytic form 
on the basis of equations of linear magneto-solid mechanics 

P u = i[5jxB], V-u = 0, (2) 

c 

5j = ^-[V x SB], SB = V x [u x B], V • B = 0. (3) 

4.7T 

Here u = u(r, t) is the field of material displacement, the fundamental dynamical variable of 
solid-mechanical theory of elastically deformable (non-flowing) material continua. These equations 
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describe Lorentz-force-d riven non-compressional vibrations of a perfectly conducting elastic mat- 
ter of a non- convective solid star with Ampere form of fluctuating current density 8j. Equation 
for SB describing coupling between fluctuating field of material displacements u and background 
magnetic field B pervading stellar material is the mathematical form of Alfven theorem about 
frozen-in lines of magnetic field in perfectly conducting matter. The adopted for Alfven vibra- 
tions of non-convective solid stars terminology [which is not new of course, see for instance[13j, 
namely, equations of magneto-solid mechanics and/or solid-magnetics is used as a solid-mechanical 
counterpart of well-known terms like equations of magneto-fluid mechanics, magnetohydrodynamics 
(MHD) and hydromagnetics pSv = (l/c)[<5j x B], <5j = (c/4vr)[V x SB], SB = V x [8v x B], 
describing flowing magento-active plasma in terms of the velocity of fluctuating flow 5v = u and 
fluctuating magnetic field SB. The MHD approach is normally utilized in astrophysics of convective 
main-sequence (MS) liquid stars such, for instance, as rapidly oscillating Ap (roAp) stars, chemically 
peculiar magnetic stars exhibiting high-frequency oscillations that have been and still remain the 
subject of extensive investigational H5] . 

Inserting (ED i n (ED the former equation of solid-magnetics takes the form 

pu{r,t) = i-[V x [V x [u(r,t) X B(r)]]] x B(r). (4) 

The analogy between perfectly conducting medium pervaded by magnetic field (magneto-active 
plasma) and elastic solid, regarded as a material continuum, is strengthened by the following ten- 
sor representation of the last equation piii = VkSM ik , where 5M ik = (l/A / K)[B i 8B k + B k 5Bi — 
Bj5Bj5ik] is the Maxwellian tensor of magnetic field stresses with SBi = Vk[uiBk — u^B^. This 
form is identical in appearance to canonical equation of solid-mechanics pxn = Vk&ik, where = 
2puik + [k — (2/3)/i] UjjSik is the Hookean tensor of mechanical stresses and = (l/2)[Vjtt& + 
VfcwJ is the tensor of shear deformations in an isotropic elastic continuous matter with shear 
modulus p and bulk modulus k (having physical dimension of pressure). The most prominent man- 
ifest of physical similarity between these two material continua is their capability of transmitting 
non-compressional perturbation by transverse waves. Unlike incompressible liquid, an elastic solid 
can respond to impulsive non-compressional load by transverse waves of shear material displace- 
ments traveling with the speed c t = yp~Jp- The unique feature of an incompressible perfectly 
conducting and magnetized continuous matter (in liquid or solid aggregated state) is the capability 
of transmitting perturbation by transverse magneto-mechanical, Alfven, wave in which mechanical 
displacements of material and fluctuations of magnetic field undergo coupled oscillations traveling 
with the speed va = B/^Arrp = a/ 2P~b / P (where Pb = B 2 /8ir is magnetic field pressure) along 
magnetic axis. It is stated, therefore, that magnetic field pervading perfectly conducting medium 
imparts to it a supplementary portion of solid-mechanical elasticity [TBI [T7]. This suggests that 
hydromagnetic Alfven vibrations of a spherical mass of a perfectly conducting matter with frozen-in 
magnetic field can be specified in a manner of eigenstates of elastic vibrations of a solid sphere. As 
for the general asteroseismology of compact objects is concerned, the above equations seems to be 
appropriate not only for neutron stars but also white dwarfs[18j and quark stars. The superdense 
material of these latter yet hypothetical compact stars is too expected to be in solid state [T91 120] . 

Equation ()4j) serves as a basis of our further analysis. The studied in above works regime of 
node-free torsion vibrations under the action of Lorentz restoring force is of some interest in that 
the rate of differentially rotational material displacements 

u(r, t) = [w(r, t) x r], w(r, t) = [V*(r)] a(t), (5) 
V 2 X (r) = 0, X (r) = ftPi{<xxO), f t (r) = A,r l (6) 
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Figure 2: The fields of material displacements in a neutron star undergoing torsional quadrupole 
(left) and octupole (right) node-free vibrations about magnetic axis. 

has one and the same form as in torsion elastic mode of node-free vibrations under the action of 
Hooke's force of mechanical shear stresses. 

Hereafter P^cosO) stands for Legendre polynomial of degree i specifying the overtone of toroidal 
a-mode. Fig. 2 shows quadrupole and octupole overtones of such vibrations. The time-dependent 
amplitude a(t) describes temporal evolution of above vibrations; the governing equation for a(t) 
is obtained form equation fllj. The prime purpose of above works was to get some insight into 
difference between spectra of discrete frequencies of toroidal a-modes in neutron star models having 
one and the same mass M and radius R, but different shapes of constant-in-time poloidal magnetic 
fields. By use of the energy method, it was found that each specific form of spatial configuration 
of static magnetic field about axis of which the neutron star matter undergoes nodeless torsional 
oscillations is uniquely reflected in the discrete frequency spectra by form of dependence of frequency 
upon overtone £ of nodeless vibration. It worth noting that first computation of discrete spectra of 
frequencies of toroidal Alfven stellar vibrations in the standing wave-regime, has been reported by 
Chandrasekhar[23j. The extensive review of other earlier computations of discrete frequency spectra 
of a-modes, uj e — uja s i — is given in well-known review of Ledoux and Walraven[14|. 

The assumption about constant in time undisturbed magnetic field means that the internal 
magnetic field pressure, Pb, the velocity va of Alfven wave in the star bulk 




(7) 



and, hence, the frequency v A = lua/2'k (where uja = va/R) and the period P A = v A of global 
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Figure 3: The basic frequency and period of global Alfven oscillations, equations (jH]), as 
functions of magnetic field intensity in the neutron star models with indicated mass and radius. 



Alfven oscillations 

b nr p 2tt [m 

remain constant in the process of vibrations whose amplitude a(t) subjects to standard equation of 
undamped harmonic oscillatorJTUJ [IT]. The allow for viscosity of stellar material leads to exponential 
damping of amplitude, but the frequency v A and, hence, the period Pa preserve one and the same 
values as in the case of non-viscous vibrations[12]. In Fig. 3 these latter quantities are plotted as 
functions of intensity B of undisturbed poloidal magnetic field in the neutron star models with 
indicated mass M and radius R. The practical usefulness of chosen logarithmic scale in this figure is 
that it shows absolute vales of v A and Pa for global vibrations of typical in mass and radius neutron 
stars. 

In what follows, we relax the assumption about constant-in-time magnetic field and examine 
the impact of its decay on the vibration energy and period. A brief analysis of such a case has 
been given in the context of magnetic white dwarfs[21j. In this paper we present a highly extensive 
consideration of this problem in the context of neutron stars with emphasis on its relevance to the 
post-quake radiation of magnetars. 
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3 Vibration-energy powered magneto-dipole emission of 
neutron star 



In the following we consider a model of a neutron star with time-dependent intensity of homogeneous 
poloidal magnetic field which can be conveniently represented in the form 



B(r,t) = fi(t)b(r), [b r = cos9,b e = -sm6,b 4 
On account of this the equation of solid-magnetics, OH), takes the form 

B\t) 



0]. 



pu(r,t) 



47T 



[V x [V x [u(r,t) x b(r,t)]]] x b(r) 



Inserting here the following separable form of fluctuating material displacements 

u(r, t) = a(r) a(t) 

we obtain 

{p&(r)}a(t) = 2P B (t) x {[V x [V x [a(r) x b(r)]]] x b(r)}a(t). 



(9) 



(10) 



(11) 



(12) 



Scalar product of ( TT2l ) with the time-independent field of instantaneous displacements a(r) followed 
by integration over the star volume leads to equation for amplitude a(t) having the form of equation 
of oscillator with depending on time spring constant 



Ma{t) + K{t)a{t) = 0, 
M=pm e , K(t)=2P B (t)k e , 



ho 



J a(r) ■ a(r) dV, a = A t V x [r r e P e (cos 9)] , 



a r 



[b(r) x [V x [V x [a(r) x b(r)]]]]dV. 



(13) 
(14) 

(15) 
(16) 



The solution of equation of non-isochronal (non-uniform in duration) and non-stationary vibrations 
with time-dependent frequency [a(t) + co 2 (t)a(t) = where uj 2 (t) = K(t)/M\ is non-trivial and 
fairly formidable task|22j. But solution of such an equation, however, is not a prime purpose of this 
work. The main subject is the impact of depletion of magnetic-field-pressure on the total energy of 
Alfven vibrations E A = (l/2)[A^d 2 + JCa 2 ] and the discrete spectrum of frequency of the toroidal 
a-mode 



U 2 e (t) = B 2 (t)Kl K 2 = 



R 2 ' 
A-npR 2 ' 



— R 2 , 



{t 2 - 1) 



2£ + 3 
2i-\ 



> 2. 



(17) 
(18) 



It is to be stated clearly from the onset that it is not our goal here to speculate about possible 
mechanisms of neutron star demagnetization and advocate conceivable laws of magnetic field decay. 
The main purpose is to gain some insight into the effect of arbitrary law of magnetic field decay in 
quaking neutron star on period of Lorentz- force-driven torsional seismic vibrations (whose quadrupole 
and octupole overtones are pictured in Fig. 2) and radiative activity of the star brought about by 
such vibrations. In the reminder of the paper we focus on the case of torsional Alfven vibrations in 
quadrupole (£ = 2) overtone. In so doing we omit index £ putting u(t) = co£ =2 (t). 
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3.1 Magnetic-field-decay induced loss of vibration energy 
The total energy stored in quake-induced Alfven seismic vibrations of the star is given by 

E , w = « + »^>, (19) 

K,{B{t)) = u 2 {B{t))M. (20) 

Perhaps most striking consequence of the magnetic-field-pressure depletion during the post-quake 
vibrational relaxation of neutron star is that it leads to the loss of vibration energy at a rate propor- 
tional to the rate of magnetic field decay 

= ^ = AW( t )B(.)gW (21) 

2 dB dt dt 

In the model under consideration, a fairly rapid decay of magnetic field during the time of post- 
quake vibrational relaxation of the star is though of as caused, to a large extent, by coupling of the 
vibrating star with material expelled by quake. In other words, escaping material removing a part of 
magnetic flux density from the star is considered to be a most plausible reason of depletion of internal 
magnetic field pressure in the star. It seems quite likely that, contrary to viscous dissipation, the loss 
of vibration energy due to decay of magnetic field must be accompanied by coherent (non-thermal) 
electromagnetic radiation. Adhering to this supposition in the next section special consideration is 
given to the conversion of the energy of Lorentz- force-driven seismic vibrations into the energy of 
magneto-dipole emission whose flux oscillates with frequency of torsional Alfven magneto-mechanical 
vibrations of the final stage solid stars, like magnetic white dwarf and neutron stars. Fig. 4 replicates 
seismic torsional Alfven vibrations of the star which are accompanied by oscillations of lines of dipolar 
magnetic field defining the beam direction of outburst X-ray emission. 

3.2 Conversion of vibration energy into power of magneto-dipole 
radiation 

The point of departure in the study of vibration-energy powered magneto-dipole emission of the star 
(whose radiation power, V, is given by Larmor's formula) is the equation 

™M = -V{tl V(t) = ±6?(t). (22) 

Consider a model of quaking neutron star whose torsional magneto-mechanical oscillations are ac- 
companied by fluctuations of total magnetic moment preserving its initial (in seismically quiescent 
state) direction: ft = /in = constant. The total magnetic dipole moment should execute oscilla- 
tions with frequency u(t) equal to that for magneto-mechanical vibrations of stellar matter which 
are described by equation for a(t). This means that 5fj,(t) and a(t) must obey equations of similar 
form, namely 

5ji(t) +cu 2 (t)5n(t) = 0, (23) 
a(t) + uj 2 (t)a(t) = 0, u 2 (t) = B 2 (t)n 2 . (24) 
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Figure 4: Schematic view of quadrupole overtone of seismic torsional Alfven vibrations of 
neutron star with homogeneous internal and dipolar external field whose lines, defining direction 
of outburst beam, oscillate with the frequency of this seismic a-mode. 



It is easy to see that equations ( 1231) and (1241) can be reconciled if 

5fi(t) = na(t). (25) 
Then, from ( !23l) . it follows 5jl = —oj 2 iia. Given this and equating 

*^ = MJc?(t)B(t)^ (26) 
dt dt 

with 

-V = -^-y^B\t)a 2 (t) (27) 

we arrive at the equation of time evolution of magnetic field 

dB(t) . 2uV . , 

which yields the following law of its decay 

B(t)= P?L= , r- 1 = 2 7 5 2 (0). (29) 

The lifetime of magnetic field r is regarded as a parameter whose value is established from below 
given relations between the period P and its time derivative P which are taken from observations. 
Knowing from observations P{t) and P(t) and estimating r one can get information about the 
magnitude of total magnetic moment and the strength of undisturbed magnetic field. 
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Figure 5: Time evolution of Alfven frequency of Lorentz-force-driven vibrations converting 
vibration energy into power of magneto-dipole radiation. 



It is worth noting that in the model of vibration-energy powered magneto-dipole emission under 
consideration, the equation of magnetic field evolution is obtained in similar fashion as equation for 
the angular velocity does in the standard model of rotation-energy powered emission of neutron 
star. As is shown in the next section, the substantial physical difference between models of rotation- 
energy and vibration-energy powered pulsating emission of neutron stars is that in the model of 
quaking neutron star vibrating in toroidal a-mode, the elongation of period of pulses is attributed 
to magnetic field decay, whereas in canonical Pacini-Gold model of radio-pulsar the lengthening of 
period of pulses is ascribed to the slow down of the neutron star rotation|39] ST 



3.3 Lengthening of vibration period 

The immediate consequence of above line of argument is the magnetic-field-decay induced length- 
ening of vibration period 



B[t) ~ ^Tfr ^ P[t) - B(ty C ~ «' 



P(t) = P(0)[l + (t/r)]^, P(0) 



C 



5(0) ! 



P(t) 



1 P(0) 



2r[l + (*/r)]V2' 
It follows that lifetime r is determined by 



(30) 
(31) 
(32) 



P 2 (0) 

P(t)P(t) = — ^ = constant. 



(33) 



This inference of the model under consideration is demonstrated in Fig. 5. The difference between 
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periods evaluated at successive moments of time t\ = and £ 2 = t is given by 



AP(f) = 
= -P(P) 



P(t) - P(0) 

g(g) 

B{t) 



(34) 



1 - 



> 0, B(t) < B(0). 



The practical usefulness of these general relations is that they can be used as a guide in search for 
fingerprints of Alfven seismic vibrations in data on oscillating emission from quaking neutron star. 



3.4 Time evolution of the vibration amplitude 

The considered model permits exact analytic solution of equation for vibration amplitude a(t) which 
is convenient to represent as 

a(t) + co 2 (t)a(t) = 0, (35) 
to 2 (t) = w(0) = u a k. (36) 

1 + t/T 

The procedure is as follows. Let us introduce new variable s = 1 + t/r. In terms of a(s), equation 
( 1351) takes the form 

sa"(s) + 2 a{s) = 0, (3 2 = w 2 (0)r 2 = const. (37) 
This equation permits exact analytic solution|24j 

a (s) = S 1/2 {d Ji(2/3s 1 / 2 ) + C 2 Y 1 {2^s 1 ' 2 )} (38) 
where Ji(2/3s 1/2 ) and Fi(2/3s 1/2 ) are Bessel functions[25] 

i r 

J 1 (z) = - / cos(>sin0-0)GZ0, (39) 
y^s) = I / sin(zsinfl - 0)d0, 2 = 2(3s 1/2 . (40) 

The arbitrary constants C\ and C 2 can be eliminated from two conditions 

a(t = 0) = o-o, Qt(t = r) = 0. (41) 

where zero-point amplitude 

, 2E A (0) 2^(0) JT(0) 

a ° = M^(oy = -rW' u{0) = ^r { ) 

is related to the average energy ^(O) stored in torsional Alfven vibrations at initial (before magnetic 
field decay) moment of time t — 0. Before magnetic field decay, the star oscillates in the harmonic 
in time regime, that is, with amplitude a = a cos(a>(0)t), so that < a 2 >= (l/2)a!jj and < a 2 >= 
(l/2)o; 2 (0)ao. The average energy of such oscillations is given by equation Ea(0) = (1/2)M < 
a 2 > +(1/2)K(0) < a 2 >= (l/2)Mw 2 (0)a§ = (l/2)AT(0)ag which relates the energy stored in 
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Figure 6: Vibration amplitude a(t) computed at fixed C and (3 and different values of rj. This 
shows that variation of this parameter is manifested in change of magnitude of |a|, but period 
elongation is not changed. 

vibrations with vibration amplitude a . As a result, the general solution of ( 1351) can be represented 
in the form 

a(t) = C[l + (t/r)] 1/2 

x{J 1 (2/3 [1 + (t/r)] 1 / 2 ) -77^(2/3 [1 + (t/r)] 1 / 2 )}, (43) 

The vibration period lengthening in the process of vibrations is illustrated in Fig. 6 and Fig. 7, where 
we plot a(t), equation fl43j), at different values of parameters (3 and r\ pointed out in the figures. 
Fig. 7 shows that 77 is the parameter regulating magnitude of vibration amplitude, the larger 77, the 
higher amplitude. However this parameter does not affect the rate of period lengthening. As it 
is clearly seen from Fig. 6 both the elongation rate of vibration period and magnitude of vibration 
amplitude are highly sensitive to parameter /3. 

All the above shows that the magnetic-field-decay induced loss of vibration energy is substantially 
different from the vibration energy dissipation caused by shear viscosity of matter resulting in heating 
of stellar material[251 [27]. As was noted, the characteristic feature of this latter mechanism of 
vibration energy conversion into the heat (i.e., into the energy of non-coherent electromagnetic 
emission responsible for the formation of photosphere of the star) is that the frequency and, hence, 
period of vibrations are the same as in the case of viscous-free vibrations[12j. However, it is no 
longer so in the case under consideration. It follows from above that depletion of magnetic field 
pressure resulting in the loss of total energy of Alfven vibrations of the star causes its vibration 
period to lengthen at a rate proportional to the rate of magnetic field decay. 



4 Numerical analysis 

To get an idea of the magnitude of characteristic parameter of vibrations providing energy supply of 
magneto-dipole radiation of a neutron star, in Table 1 we present results of numerical computations 
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Figure 7: Vibration amplitude a(t) computed at fixed C and 77 and different values of (3. 



of the fundamental frequency of neutron star oscillations in quadrupole toroidal a-mode and time 
of decay of magnetic field r as functions of increasing magnetic field. As was emphasized, the most 
striking feature of considered model of vibration powered radiation is the lengthening of periods 
of pulsating emission caused by decay of internal magnetic field. This suggests that this model is 
relevant to electromagnetic activity of magnetars - neutron stars endowed with magnetic field of 
extremely high intensity the radiative activity of which is ultimately related to the magnetic field 
decay. Such a view is substantiated by estimates of Alfven frequency presented in the table. For 
magnetic fields of typical rotation powered radio pulsars, B ~ 10 12 G, the computed frequency va is 
much smaller than the detected frequency of pulses whose origin is attributed to lighthouse effect. 
In the meantime, for neutron stars with magnetic fields B ~ 10 14 G the estimates of 2/4 are in the 
realm of observed frequencies of high-energy pulsating emission of soft gamma repeaters (SGRs), 
anomalous X-ray pulsars (AXPs) and sources exhibiting similar features. According to common 
belief, these are magnetars - highly magnetized neutron stars whose radiative activity is related with 
magnetic field decay. The amplitude of vibrations is estimated as 

1/2 

= 3.423 x lO^E^B^R- 3/2 . (45) 

where £^,40 = -Ea/(10 40 erg) is the energy stored in the vibrations. R 6 = R/(10 6 cm) and 
B14 = £?/(10 14 G). The presented computations show that the decay time of magnetic field (equal 
to duration time of vibration powered radiation in question) strongly depends on the intensity of 
initial magnetic field of the star: the larger magnetic field B the shorter time of radiation r at the 
expense of energy of vibration in decay during this time magnetic field. The effect of equation of 
state of neutron star matter (which is most strongly manifested in different values mass and radius 
of the star) on frequency va is demonstrated by numerical vales of this quantity for magnetars with 
one and the same value of magnetic field B = 10 15 G but different values of mass and radius. 



2E A (0) 
Moo 2 (0) 
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Table 1: The Alfven frequency of Lorentz- force-driven torsion vibrations, va, arid their lifetime 
equal to decay time of magnetic field, r, in neutron stars with magnetic fields typical for pulsars 
and magnet ars. 





M(M ) 


R{km) 


B(G) 


z/ A (Hz) 


7"(yr) 


Pulsars 


0.8 


20 


10 12 


3.25 x 10" 3 


4.53 x 10 10 




1.0 


15 


10 13 


2.52 x 10~ 2 


2.98 x 10 7 


Magnet ars 


1.1 


13 


10 14 


0.22 


7.4 x 10 3 




1.2 


12 


10 15 


2.06 


1.31 




1.3 


11 


10 15 


1.89 


2.38 




1.4 


10 


10 16 


17.4 


4.44 x 10~ 4 



The luminosity powered by neutron star vibrations in quadrupole toroidal a-mode is given by 

V = ^B\t)a 2 {t), fx = {1/2)B{0)R 3 , B(t) = B(0)[1 + t/r}- 1 ' 2 , (46) 
a(t) = Cs 1/2 {J 1 (z(t))-r ] Y 1 (z(t))} J z = 2a;(0)r (1 + t/r) (47) 

The presented in Fig. 8 computations of power of magneto-dipole radiation of a neutron stars (of one 
and the same mass and radius but different values of magnetic fields) exhibit oscillating character 
of luminosity. The frequency of these oscillations equal to that of torsional Alfven seismic vibrations 
of neutron star. 

All above suggests that developed theory of vibration-energy powered emission of neutron star 
is relevant to electromagnetic activity of magnetars - neutron stars endowed with magnetic field 
of extremely high intensity the radiative activity of which is ultimately related to the magnetic 
field decay. This subclass of highly magnetized compact objects is commonly associated with soft 
gamma repeaters and anomalous X-ray pulsars[29l [30l [311 [32] - young isolated and seismically active 
neutron stars[33l [34]- The magnetar quakes are exhibited by short-duration thermonuclear gamma- 
ray flash followed by rapidly oscillating X-ray flare of several-hundred-seconds duration. During 
this latter stage of quake-induced radiation of magnetar, a long-periodic (2-12 sec) modulation of 
brightness was observed lasting for about 3-4 minutes. Such long period of pulsed emission might 
be expected from old rotation-energy powered pulsars (due to slow down of rotation), but not 
from magnetars which are young neutron stars, as follows from their association with pretty young 
supernovae. Taking this into account (and also the fact that energy release during X-ray flare is 
much larger than the energy of rigid-body rotation with so long periods) it has been suggested[2j 
that detected long-periodic pulsating emission of magnetars is powered by the energy of torsional 
magneto-elastic vibrations triggered by quake. In view of key role of ultra-strong magnetic field 
it is quite likely that quasi-periodic oscillations (QPOs) of outburst flux from SGR 1806-20 and 
SGR 1900+14 detected i n [351 ESI EZ] are produced by torsional seismic vibrations predominately 
sustained by Lorentz force [101 HH Q2] which are accompanied, as was argued above, by monotonic 
decay of background magnetic field. If so, the predicted elongation of QPOs period of oscillating 
outburst emission from quaking magnetars should be traced in existing and future observations. 
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Figure 8: Time evolution of luminosity of magneto-dipole radiation powered by energy of 
torsional Alfven seismic vibrations of neutron star with mass M = 1.2M and radius R = 15 
km with intensities of magnetic field. 

4.1 Comparison with the rotation powered neutron star 

For a sake of comparison, in the considered model of vibration powered radiation, the equation of 
magnetic field evolution is obtained in similar fashion as that for the angular velocity Q(t) does in 
the standard model of rotation powered neutron star which rests on equations [38J 

E R (t) = \nnt), I = \mW (49) 
2 5 

One of the basic postulates of the model of rotation-energy powered emission is that the time 
evolution of total magnetic moment of the star is governed by the equation 

Sfl(t) = [ft(t) x [n(i) x //]], /x = constant. (50) 

It follows 

5jx 2 (t) = /ii^ 4 (t), fi± = \i sin# (51) 

where 8 is angle of inclination of fi to O(t). The total magnetic moment of non-rotating neutron 
star is parametrized by equation of uniformly magnetized (along the polar axis) sphere 

H = /in, fi = ^BR Z = constant. (52) 

This parametrization presumes that in the rotation powered neutron star, the frozen-in the star 
magnetic field operates like a passive promoter of magneto- dipole radiation, that is, intensity of 
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the internal magnetic field remain constant in the process of radiation. As a result, the equation of 
energy conversion from rotation to magnetic dipole radiation is reduced to equation of slow down 
of rotation 

n(t) = -K{l 3 (t), K = ^ t (53) 
fi(t) = -£E=, T-^ = 2Ktf(0). (54) 

V 1 + t / T 

where 9 is angle of inclination of \x to Ct(t). From the last equation it follows 

F(«) = F(0)[l + (i/T)] 1/2 , P(0) = ^-y (55) 

= iw^MF (56) 

and, hence, the lifetime r is related with P{t) and P(t) as 

P 2 (0) 

P(t)P(t) = — — = constant. (57) 
2r 

Equating two independent estimates for r given by equations ( l54l) and ( l57f ) we arrive at widely 
utilized analytic estimate of magnetic field on the neutron star pole: 



B = [3/c 3 /(27r 2 i? 6 )] 1/2 y / P(t)P(t). 
For a neutron star of mass M = M , and radius R = 13 km, one has 



B = 3.210 19 ^/P(t)P(t), G. 

The outlined treatment of rotation-powered magneto-dipole radiation of a neutron star empha- 
sizes kinematic nature of variation of magnetic moment of pulsar whose magnetic field is regarded 
as independent of time. Thus, the substantial physical difference between vibration-powered and 
rotation-powered neutron star models is that in the former the elongation of pulse period is at- 
tributed to magnetic field decay, whereas in the latter the period lengthening is ascribed to slow 
down of rotation [391 1401 l4Tj . 

4.2 Comment on magnetic field decay in quaking neutron star 

The considered law of magnetic field decay cannot be, of course, regarded as universal because it 
reflects a quite concrete line of argument regarding the fluctuations of magnetic moment of the 
star. The interrelation between quake-induced oscillations of total magnetic moment of the star, 
5fj,(t), and the amplitude, a(t), of its seismic magneto-mechanical oscillations can be consistently 
interpreted with the aid of the function of dipole demagnetization f(_B(t)) which is defined by the 
following condition of self-consistency in a of right and left hand sides of equation ( 1221) . namely 

5jl(t) = f(B(t))a(t). (58) 
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The vector-function of dipole demagnetization f(B(t)) depending on decaying magnetic field re- 
flects temporal changes of electromagnetic properties of neutron star matter as well as evolution 
of magnetic-field-promoted coupling between neutron star and its environment. This means that 
specific form of this phenomenological function should be motivated by heuristic arguments taking 
into account these factors. With this form of 5jl(t), equation ((22D is transformed to magnetic field 
decay of the form 



dB(t) 



-v 



f\B(t)) 



V 



const. 



dt '' B(t) ' ' 3Mk 2 c 3 

In the above considered case, this function is given by 

f (B(t)) = B{t) B(t), = k 2 /i = constant. 



(59) 



(60) 



The practical usefulness of the dipole demagnetization function, f(B(t)), consists in that it provides 
economic way of studying a vast variety of heuristically motivated laws of magnetic field decay, 
B = B{t), whose inferences can ultimately be tested by observations. In this subsection, with no 
discussing any specific physical mechanism which could be responsible for magnetic field decay, we 
consider a set of representative examples of demagnetization function f (£(£)) some of which have 
been regarded before, though in a somewhat different context [42-49]. Here we stress again that the 
model under consideration deals with magnetic field decay in the course of vibrations triggered by 
starquake, not with long-term secular decay which has been the subject of these latter investigations 
(see also references therein). 
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: P(l)=P(0)exp((/r) 

:" ' ' ' ' I i ' ' ' I ' ' 


example 1 


f(i(0)=A/»p(i) 


' ' i ' i i i i ' i ' ' 


[ P(0 = -P(0)(l+;/r) 






example 2 


— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1- 

Vr + I y,B(l) 

: P(t)=P(0)[l-ln(l+(/T)] 4 


' ' i ' i i i i ' i ' ' 
example 3 



" a 

2.0 S- 



0.0 0.2 0.4 0.6 0.8 



0.2 0.4 OS 0.8 



Figure 9: Three representative examples with different laws of the magnetic field decay resulting 
in the elongation of the vibration period for each case. 

1. As a first representative example, a model of quaking neutron star whose function of dipole 
demagnetization has the following form 

f(B(f)) = 0B(t) (61) 
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is been considered. In such demagnetization model, the temporal evolution of the magnetic field 
and the lengthening of the vibration period obey the following laws 



B(t) = B(0)e~ t/T , P(t) = P(0)e t/r , r 



m 
p( t y 



(62) 



where B(0) is the intensity of magnetic field before quake. 2. For a quaking neutron star model 
whose function of demagnetization is given by 

f(B(t)) = Py/B®B(t) (63) 

the resultant equation of magnetic field evolution and the vibration period elongation read 

-l 



B(t) = P(0) 1 + 



t 



(64) 



Similar analysis can be performed for the demagnetization function of the form 
Namely, 



f (£(*)) = /3y/B™- 1 {t)B(t), m = 3,4,5... 
dB(t) 



dt 
B(t) = 



= - 1 B m (t), 1 = v f3 2 , 
B(0) 



TTT' r-^m-l)^- 1 ^). 



[l+t/Tm] l/(m-l)' 

3. Finally, let's consider a model with quite sophisticated function of dipole demagnetization 



f(B(t)) 



B(0) B(t) 

VvW) 



T + t 



(65) 



which lead to a fairly non-trivial logarithmic law of magnetic field decay and vibration period length- 
ening 



Bit) = B(0) 
P(t) = P(0) 



1 - In 1 + 



1 - In ( 1 + - 



~i -i 



T = 



P(0) 
P(0)' 



(66) 
(67) 



For each of above examples, the magnetic field and the resultant lengthening of vibration period, 
computed as functions of fractional time t/r, are shown in Fig. 9. This third example is interesting in 
that computed in this model ratio z>(t)/z/(0), pictured in Fig. 10, as a function of t, is similar to that 
which exhibit data on post-glitch emission of PSR J1846-0258 [50]. These examples show that the 
period elongation is the common effect of magnetic field decay. The interrelations between periods 
and its derivatives substantially depend on specific form of the magnetic field decay law although 
physical processes responsible for magnetic field decay remain uncertain. Practical significance of 
above heuristic line of argument is that it leads to meaningful conclusion (regarding elongation of 
periods of oscillating magneto-dipole emission) even when detailed mechanisms of magnetic field 
decay in the course of vibrations are not exactly known. 
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Figure 10: The rate of frequency z>(x) normalized to u(0) as a function of x = t/r for 
logarithmic law of magnetic filed decay. 

5 Summary 

It is generally realized today that the standard model of inclined rotator, lying at the base of our 
understanding of radio pulsars, faces serious difficulties in explaining the long-periodic (2 < P < 
12 s) pulsed radiation of soft gamma repeaters (SGRs) and anamalous X-ray pulsars (AXPs). 
Observations show that persistent X-ray luminosity of these sources (10 34 < L x < 10 36 erg s _1 ) 
is appreciably (10-100 times) larger than expected from neutron star deriving radiation power from 
energy of rotation with frequency of detected pulses. It is believed that this discrepancy can be 
resolved assuming that AXP/SGR-like sources are magnetars - young, isolated and seismically active 
neutron stars whose energy supply of pulsating high-energy radiation comes not from rotation (as is 
the case of radio pulsars) but from different process involving decay of ultra strong magnetic field, 
10 14 < B < 10 16 G. Adhering to this attitude we have presented the model of quaking neutron 
star deriving radiation power from the energy of torsional Lorentz-force-d riven oscillations. It is 
appropriate to remind early works of the infancy of neutron star era |5TJ [55] in which it has been 
pointed out for the first time that vibrating neutron star should operate like Hertzian magnetic 
dipole deriving radiative power of magneto-dipole emission from the energy of magneto-mechanical 
vibrations[53]. What is newly disclosed here is that the main prerequisite of the energy conversion 
from vibrations into radiation is the decay of magnetic field in the star. Since the magnetic field 
decay is one of the most conspicuous features distinguishing magnetars from rotation powered 
pulsars, it seems meaningful to expect that at least some of AXP/SGR - like sources are magnetars 
deriving power of pulsating magnetic dipole radiation from the energy of quake-induced torsion 
seismic vibrations. 
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